Analysis of intra-specific variations in the venom of individual snakes based on Raman spectroscopy.

The knowledge of variations in the composition of venoms from different snakes is important from both theoretical and practical points of view, in particular, at developing and selecting an antivenom. Many studies on this topic are conducted with pooled venoms, while the existence and significance of variations in the composition of venoms between individual snakes of the same species are emphasized by many authors. It is important to study both inter- and intra-specific, including intra-population, venom variations, because intra-specific variations in the venom composition may affect the effectiveness of antivenoms as strongly as inter-specific. In this work, based on venom Raman spectroscopy with principal component analysis, we assessed the variations in venoms of individual snakes of the Vipera nikolskii species from two populations and compared these intra-specific variations with inter-specific variations (with regard to the other related species). We demonstrated intra-specific (inter- and intra-population) differences in venom compositions which are smaller than inter-specific variations. We also assessed the compositions of V. nikolskii venoms from two populations to explain inter-population differences. The method used is rapid and requires virtually no preparation of samples, used in extremely small quantities, allowing the venoms of individual snakes to be analyzed. In addition, the method is informative and capable of detecting fairly subtle differences in the composition of venoms.

Effects of the Heterodimeric Neurotoxic Phospholipase A2 from the Venom of Vipera nikolskii on the Contractility of Rat Papillary Muscles and Thoracic Aortas.

Phospholipases A2 (PLA2s) are a large family of snake toxins manifesting diverse biological effects, which are not always related to phospholipolytic activity. Snake venom PLA2s (svPLA2s) are extracellular proteins with a molecular mass of 13-14 kDa. They are present in venoms in the form of monomers, dimers, and larger oligomers. The cardiovascular system is one of the multiple svPLA2 targets in prey organisms. The results obtained previously on the cardiovascular effects of monomeric svPLA2s were inconsistent, while the data on the dimeric svPLA2 crotoxin from the rattlesnake Crotalus durissus terrificus showed that it significantly reduced the contractile force of guinea pig hearts. Here, we studied the effects of the heterodimeric svPLA2 HDP-1 from the viper Vipera nikolskii on papillary muscle (PM) contractility and the tension of the aortic rings (ARs). HDP-1 is structurally different from crotoxin, and over a wide range of concentrations, it produced a long-term, stable, positive inotropic effect in PMs, which did not turn into contractures at the concentrations studied. This also distinguishes HDP-1 from the monomeric svPLA2s, which at high concentrations inhibited cardiac function. HDP-1, when acting on ARs preconstricted with 10 µM phenylephrine, induced a vasorelaxant effect, similar to some other svPLA2s. These are the first indications of the cardiac and vascular effects of true vipers' heterodimeric svPLA2s.

Differentiation of snake venom using Raman spectroscopic analysis.

Snake venoms are complex mixtures of different substances, proteins being their predominant components. To study the composition of venoms, methods based on chromatographic separation and mass spectrometric analysis are currently used, requiring the application of a number of sophisticated instruments. To assess the composition of snake venoms, we propose an alternative method based on Raman spectroscopy, which is an express method to study the structural features of different substances, including proteins. The method does not require preliminary preparation of the samples, which are used in small quantities; this makes Raman spectroscopy extremely attractive for venom research. In this work, we have carried out Raman spectroscopic studies on a number of dry venoms from various venomous snakes. Based on the obtained Raman spectra, with the help of mathematical methods of dimensionality reduction and clustering, differentiation of venoms reflecting their composition and the assignment of the venom producing snake to the corresponding family or even genus were performed. The proposed method can be used to analyze both the composition of and variations in venoms of different snake species, including rare and endangered ones.

Effects of Cobra Cardiotoxins on Intracellular Calcium and the Contracture of Rat Cardiomyocytes Depend on Their Structural Types.

Cardiotoxins (CaTx) of the three-finger toxin family are one of the main components of cobra venoms. Depending on the structure of the N-terminal or the central polypeptide loop, they are classified into either group I and II or P- and S-types, respectively, and toxins of different groups or types interact with lipid membranes variably. While their main target in the organism is the cardiovascular system, there is no data on the effects of CaTxs from different groups or types on cardiomyocytes. To evaluate these effects, a fluorescence measurement of intracellular Ca2+ concentration and an assessment of the rat cardiomyocytes' shape were used. The obtained results showed that CaTxs of group I containing two adjacent proline residues in the N-terminal loop were less toxic to cardiomyocytes than group II toxins and that CaTxs of S-type were less active than P-type ones. The highest activity was observed for Naja oxiana cobra cardiotoxin 2, which is of P-type and belongs to group II. For the first time, the effects of CaTxs of different groups and types on the cardiomyocytes were studied, and the data obtained showed that the CaTx toxicity to cardiomyocytes depends on the structures both of the N-terminal and central polypeptide loops.

The Potassium Channel Blocker ß-Bungarotoxin from the Krait Bungarus multicinctus Venom Manifests Antiprotozoal Activity.

Protozoal infections are a world-wide problem. The toxicity and somewhat low effectiveness of the existing drugs require the search for new ways of protozoa suppression. Snake venom contains structurally diverse components manifesting antiprotozoal activity; for example, those in cobra venom are cytotoxins. In this work, we aimed to characterize a novel antiprotozoal component(s) in the Bungarus multicinctus krait venom using the ciliate Tetrahymena pyriformis as a model organism. To determine the toxicity of the substances under study, surviving ciliates were registered automatically by an original BioLaT-3.2 instrument. The krait venom was separated by three-step liquid chromatography and the toxicity of the obtained fractions against T. pyriformis was analyzed. As a result, 21 kDa protein toxic to Tetrahymena was isolated and its amino acid sequence was determined by MALDI TOF MS and high-resolution mass spectrometry. It was found that antiprotozoal activity was manifested by ß-bungarotoxin (ß-Bgt) differing from the known toxins by two amino acid residues. Inactivation of ß-Bgt phospholipolytic activity with p-bromophenacyl bromide did not change its antiprotozoal activity. Thus, this is the first demonstration of the antiprotozoal activity of ß-Bgt, which is shown to be independent of its phospholipolytic activity.

S- and P-type cobra venom cardiotoxins differ in their action on isolated rat heart.

Background: The cardiovascular system is one of the first systems to be affected by snake toxins; but not many toxins exert a direct effect on the heart. Cobra venom cardiotoxins are among those few toxins that attack the heart. Although the two cardiotoxin types (S and P) differ in their central-loop structure, it is not known whether they differ in their effect on the mammalian heart. We compared the effects of S- and P-type cardiotoxins, CTÐ¥-1 and CTÐ¥-2, respectively, from the cobra Naja oxiana, on the isolated rat heart. Methods: An isolated rat heart perfused according to the Langendorff technique was used in this study to investigate the activity of cardiotoxins CTX-1 and CTX-2. The following parameters were registered: the left ventricular developed pressure, calculated as the difference between systolic and diastolic pressure in the left ventricle, the end-diastolic pressure, the heart rate, time to maximal end-diastolic pressure (heart contracture), and time to depression of the heart contraction. Results: Both cardiotoxins at the concentration of 5 µg/mL initially produce a slight increase in systolic intraventricular pressure, followed by its rapid decrease with a simultaneous increase in diastolic intraventricular pressure until reaching contracture. CTX-2 blocks cardiac contractions faster than CTX-1; in its presence the maximum diastolic pressure is reached faster and the magnitude of the developed contracture is higher. Conclusion: The P-type cardiotoxin CTX-2 more strongly impairs rat heart functional activity than the S-type cardiotoxin CTX-1, as expressed in its faster blockage of cardiac contractions as well as in more rapid development and greater magnitude of contracture in its presence.

Anti-HIV Activity of Snake Venom Phospholipase A2s: Updates for New Enzymes and Different Virus Strains.

Since the beginning of the HIV epidemic, lasting more than 30 years, the main goal of scientists was to develop effective methods for the prevention and treatment of HIV infection. Modern medicines have reduced the death rate from AIDS by 80%. However, they still have side effects and are very expensive, dictating the need to search for new drugs. Earlier, it was shown that phospholipases A2 (PLA2s) from bee and snake venoms block HIV replication, the effect being independent on catalytic PLA2 activity. However, the antiviral activity of human PLA2s against Lentiviruses depended on catalytic function and was mediated through the destruction of the viral membrane. To clarify the role of phospholipolytic activity in antiviral effects, we analyzed the anti-HIV activity of several snake PLA2s and found that the mechanisms of their antiviral activity were similar to that of mammalian PLA2. Our results indicate that snake PLA2s are capable of inhibiting syncytium formation between chronically HIV-infected cells and healthy CD4-positive cells and block HIV binding to cells. However, only dimeric PLA2s had pronounced virucidal and anti-HIV activity, which depended on their catalytic activity. The ability of snake PLA2s to inactivate the virus may provide an additional barrier to HIV infection. Thus, snake PLA2s might be considered as candidates for lead molecules in anti-HIV drug development.

Variability in the Spatial Structure of the Central Loop in Cobra Cytotoxins Revealed by X-ray Analysis and Molecular Modeling.

Cobra cytotoxins (CTs) belong to the three-fingered protein family and possess membrane activity. Here, we studied cytotoxin 13 from Naja naja cobra venom (CT13Nn). For the first time, a spatial model of CT13Nn with both "water" and "membrane" conformations of the central loop (loop-2) were determined by X-ray crystallography. The "water" conformation of the loop was frequently observed. It was similar to the structure of loop-2 of numerous CTs, determined by either NMR spectroscopy in aqueous solution, or the X-ray method. The "membrane" conformation is rare one and, to date has only been observed by NMR for a single cytotoxin 1 from N. oxiana (CT1No) in detergent micelle. Both CT13Nn and CT1No are S-type CTs. Membrane-binding of these CTs probably involves an additional step-the conformational transformation of the loop-2. To confirm this suggestion, we conducted molecular dynamics simulations of both CT1No and CT13Nn in the Highly Mimetic Membrane Model of palmitoiloleoylphosphatidylglycerol, starting with their "water" NMR models. We found that the both toxins transform their "water" conformation of loop-2 into the "membrane" one during the insertion process. This supports the hypothesis that the S-type CTs, unlike their P-type counterparts, require conformational adaptation of loop-2 during interaction with lipid membranes.

Effects of Cardiotoxins from Naja oxiana Cobra Venom on Rat Heart Muscle and Aorta: A Comparative Study of Toxin-Induced Contraction Mechanisms.

Cardiotoxins (CaTxs) are a group of snake toxins that affect the cardiovascular system (CVS). Two types (S and P) of CaTxs are known, but the exact differences in the effects of these types on CVS have not been thoroughly studied. We investigated cellular mechanisms of action on CVS for Naja oxiana cobra CaTxs CTX-1 (S-type) and CTX-2 (P-type) focusing on the papillary muscle (PM) contractility and contraction of aortic rings (AR) supplemented by pharmacological analysis. It was found that CTX-1 and CTX-2 exerted dose-dependent effects manifested in PM contracture and AR contraction. CTX-2 impaired functions of PM and AR more strongly than CTX-1. Effects of CaTxs on PM were significantly reduced by nifedipine, an L-type Ca2+ channel blocker, and by KB-R7943, an inhibitor of reverse-mode Na+/Ca2+ exchange. Furthermore, 2-aminoethoxydiphenyl borate, an inhibitor of store-operated calcium entry, partially restored PM contractility damaged by CaTxs. The CaTx influence on AR contracture was significantly reduced by nifedipine and KB-R7943. The involvement of reverse-mode Na+/Ca2+ exchange in the effect of CaTxs on the rat aorta was shown for the first time. The results obtained indicate that CaTx effects on CVS are mainly associated with disturbance of transporting systems responsible for the Ca2+ influx.

Membrane-Disrupting Activity of Cobra Cytotoxins Is Determined by Configuration of the N-Terminal Loop.

In aqueous solutions, cobra cytotoxins (CTX), three-finger folded proteins, exhibit conformational equilibrium between conformers with either cis or trans peptide bonds in the N-terminal loop (loop-I). The equilibrium is shifted to the cis form in toxins with a pair of adjacent Pro residues in this loop. It is known that CTX with a single Pro residue in loop-I and a cis peptide bond do not interact with lipid membranes. Thus, if a cis peptide bond is present in loop-I, as in a Pro-Pro containing CTX, this should weaken its lipid interactions and likely cytotoxic activities. To test this, we have isolated seven CTX from Naja naja and N. haje cobra venoms. Antibacterial and cytotoxic activities of these CTX, as well as their capability to induce calcein leakage from phospholipid liposomes, were evaluated. We have found that CTX with a Pro-Pro peptide bond indeed exhibit attenuated membrane-perturbing activity in model membranes and lower cytotoxic/antibacterial activity compared to their counterparts with a single Pro residue in loop-I.

Snake venom phospholipase A2s exhibit strong virucidal activity against SARS-CoV-2 and inhibit the viral spike glycoprotein interaction with ACE2.

The COVID-19 pandemic caused by SARS-CoV-2 requires new treatments both to alleviate the symptoms and to prevent the spread of this disease. Previous studies demonstrated good antiviral and virucidal activity of phospholipase A2s (PLA2s) from snake venoms against viruses from different families but there was no data for coronaviruses. Here we show that PLA2s from snake venoms protect Vero E6 cells against SARS-CoV-2 cytopathic effects. PLA2s showed low cytotoxicity to Vero E6 cells with some activity at micromolar concentrations, but strong antiviral activity at nanomolar concentrations. Dimeric PLA2 from the viper Vipera nikolskii and its subunits manifested especially potent virucidal effects, which were related to their phospholipolytic activity, and inhibited cell-cell fusion mediated by the SARS-CoV-2 spike glycoprotein. Moreover, PLA2s interfered with binding both of an antibody against ACE2 and of the receptor-binding domain of the glycoprotein S to 293T/ACE2 cells. This is the first demonstration of a detrimental effect of PLA2s on ß-coronaviruses. Thus, snake PLA2s are promising for the development of antiviral drugs that target the viral envelope, and could also prove to be useful tools to study the interaction of viruses with host cells.

The omega-loop of cobra cytotoxins tolerates multiple amino acid substitutions.

Cobra cytotoxins (CTs), the three-fingered proteins, feature high amino acid sequence homology in the beta-strands and variations in the loop regions. We selected a pair of cytotoxins from Naja kaouthia crude venom to clarify the sequence-structure relationships. Using chromatography and mass spectroscopy, we separated and identified the mixture of cytotoxins 2 and 3, differentiated by the only Val 41/Ala 41 substitution. Here, using natural abundance 13C, 15N NMR-spectroscopy we performed chemical shift assignments of the signals of the both toxins in aqueous solution in the major and minor forms. Combining NOE and chemical shift data, the toxins' spatial structure was determined. Finally, we proved that the tip of the "finger"-2, or the loop-2 of cytotoxins adopts the shape of an omega-loop with a tightly-bound water molecule in its cavity. Comparison with other NMR and X-ray structures of cytotoxins possessing different amino acid sequences reveals spatial similarity in this family of proteins, including the loop-2 region, previously considered to be flexible.

Snake C-type lectin-like proteins inhibit nicotinic acetylcholine receptors.

Venoms of viperid snakes affect mostly hemostasis, while C-type lectin-like proteins (CTLPs), one of the main components of viperid venoms, act as anticoagulants, procoagulants, or agonists/antagonists of platelet activation. However, we have shown earlier that CTLPs from the saw-scaled viper Echis multisquamatus, called emunarecins EM1 and EM2, were able to inhibit nicotinic acetylcholine receptors (nAChRs) in neurons of a pond snail (Lymnaea stagnalis). Here we analysed the structure of the emunarecins by mass spectrometry and report that EM1 and EM2 inhibit fluorescent α-bungarotoxin binding to both muscle-type nAChRs from Torpedo californica and human neuronal α7 nAChRs. EM1 at 23µM and EM2 at 9µM almost completely prevented fluorecsent α-bungarotoxin binding to muscle-type nAChRs. Interaction with human neuronal α7 nAChR was weaker; EM1 at the concentration of 23µM blocked the α-bungarotoxin binding only by about 40% and EM2 at 9µM by about 20%. The efficiency of the EM2 interaction with nAChRs was comparable to that of a non-conventional toxin, WTX, from Naja kaouthia cobra venom. Together with the data obtained earlier, these results show that CTLPs may represent new nAChR ligands.

Novel Bradykinin-Potentiating Peptides and Three-Finger Toxins from Viper Venom: Combined NGS Venom Gland Transcriptomics and Quantitative Venom Proteomics of the Azemiops feae Viper.

Feae's viper Azemipos feae belongs to the Azemiopinae subfamily of the Viperidae family. The effects of Viperidae venoms are mostly coagulopathic with limited neurotoxicity manifested by phospholipases A2. From A. feae venom, we have earlier isolated azemiopsin, a novel neurotoxin inhibiting the nicotinic acetylcholine receptor. To characterize other A. feae toxins, we applied label-free quantitative proteomics, which revealed 120 unique proteins, the most abundant being serine proteinases and phospholipases A2. In total, toxins representing 14 families were identified, among which bradykinin-potentiating peptides with unique amino acid sequences possessed biological activity in vivo. The proteomic analysis revealed also basal (commonly known as non-conventional) three-finger toxins belonging to the group of those possessing neurotoxic activity. This is the first indication of the presence of three-finger neurotoxins in viper venom. In parallel, the transcriptomic analysis of venom gland performed by Illumina next-generation sequencing further revealed 206 putative venom transcripts. Together, the study unveiled the venom proteome and venom gland transciptome of A. feae, which in general resemble those of other snakes from the Viperidae family. However, new toxins not found earlier in viper venom and including three-finger toxins and unusual bradykinin-potentiating peptides were discovered.

Screening Snake Venoms for Toxicity to Tetrahymena Pyriformis Revealed Anti-Protozoan Activity of Cobra Cytotoxins.

Snake venoms possess lethal activities against different organisms, ranging from bacteria to higher vertebrates. Several venoms were shown to be active against protozoa, however, data about the anti-protozoan activity of cobra and viper venoms are very scarce. We tested the effects of venoms from several snake species on the ciliate Tetrahymena pyriformis. The venoms tested induced T. pyriformis immobilization, followed by death, the most pronounced effect being observed for cobra Naja sumatrana venom. The active polypeptides were isolated from this venom by a combination of gel-filtration, ion exchange and reversed-phase HPLC and analyzed by mass spectrometry. It was found that these were cytotoxins of the three-finger toxin family. The cytotoxins from several cobra species were tested and manifested toxicity for infusorians. Light microscopy revealed that, because of the cytotoxin action, the infusorians' morphology was changed greatly, from teardrop-like to an almost spherical shape, this alteration being accompanied by a leakage of cell contents. Fluorescence microscopy showed that the fluorescently labelled cytotoxin 2 from cobra N. oxiana was localized mainly at the membrane of killed infusorians, indicating that cytotoxins may kill T. pyriformis by causing membrane rupture. This work is the first evidence of the antiprotozoal activity of cobra venom cytotoxins, as demonstrated by the example of the ciliate T. pyriformis.

Pancreatic and snake venom presynaptically active phospholipases A2 inhibit nicotinic acetylcholine receptors.

Phospholipases A2 (PLA2s) are enzymes found throughout the animal kingdom. They hydrolyze phospholipids in the sn-2 position producing lysophospholipids and unsaturated fatty acids, agents that can damage membranes. PLA2s from snake venoms have numerous toxic effects, not all of which can be explained by phospholipid hydrolysis, and each enzyme has a specific effect. We have earlier demonstrated the capability of several snake venom PLA2s with different enzymatic, cytotoxic, anticoagulant and antiproliferative properties, to decrease acetylcholine-induced currents in Lymnaea stagnalis neurons, and to compete with α-bungarotoxin for binding to nicotinic acetylcholine receptors (nAChRs) and acetylcholine binding protein. Since nAChRs are implicated in postsynaptic and presynaptic activities, in this work we probe those PLA2s known to have strong presynaptic effects, namely ß-bungarotoxin from Bungarus multicinctus and crotoxin from Crotalus durissus terrificus. We also wished to explore whether mammalian PLA2s interact with nAChRs, and have examined non-toxic PLA2 from porcine pancreas. It was found that porcine pancreatic PLA2 and presynaptic ß-bungarotoxin blocked currents mediated by nAChRs in Lymnaea neurons with IC50s of 2.5 and 4.8 µM, respectively. Crotoxin competed with radioactive α-bungarotoxin for binding to Torpedo and human α7 nAChRs and to the acetylcholine binding protein. Pancreatic PLA2 interacted similarly with these targets; moreover, it inhibited radioactive α-bungarotoxin binding to the water-soluble extracellular domain of human α9 nAChR, and blocked acetylcholine induced currents in human α9α10 nAChRs heterologously expressed in Xenopus oocytes. These and our earlier results show that all snake PLA2s, including presynaptically active crotoxin and ß-bungarotoxin, as well as mammalian pancreatic PLA2, interact with nAChRs. The data obtained suggest that this interaction may be a general property of all PLA2s, which should be proved by further experiments.

Cobra Venom Factor and Ketoprofen Abolish the Antitumor Effect of Nerve Growth Factor from Cobra Venom.

We showed recently that nerve growth factor (NGF) from cobra venom inhibited the growth of Ehrlich ascites carcinoma (EAC) inoculated subcutaneously in mice. Here, we studied the influence of anti-complementary cobra venom factor (CVF) and the non-steroidal anti-inflammatory drug ketoprofen on the antitumor NGF effect, as well as on NGF-induced changes in EAC histological patterns, the activity of lactate and succinate dehydrogenases in tumor cells and the serum level of some cytokines. NGF, CVF and ketoprofen reduced the tumor volume by approximately 72%, 68% and 30%, respectively. The antitumor effect of NGF was accompanied by an increase in the lymphocytic infiltration of the tumor tissue, the level of interleukin 1β and tumor necrosis factor α in the serum, as well as the activity of lactate and succinate dehydrogenases in tumor cells. Simultaneous administration of NGF with either CVF or ketoprofen abolished the antitumor effect and reduced all other effects of NGF, whereas NGF itself significantly decreased the antitumor action of both CVF and ketoprofen. Thus, the antitumor effect of NGF critically depended on the status of the immune system and was abolished by the disturbance of the complement system; the disturbance of the inflammatory response canceled the antitumor effect as well.

Peptides from puff adder Bitis arietans venom, novel inhibitors of nicotinic acetylcholine receptors.

Phospholipase A2 (named bitanarin) possessing capability to block nicotinic acetylcholine receptors (nAChRs) was isolated earlier (Vulfius et al., 2011) from puff adder Bitis arietans venom. Further studies indicated that low molecular weight fractions of puff adder venom inhibit nAChRs as well. In this paper, we report on isolation from this venom and characterization of three novel peptides called baptides 1, 2 and 3 that reversibly block nAChRs. To isolate the peptides, the venom of B. arietans was fractionated by gel-filtration and reversed phase chromatography. The amino acid sequences of peptides were established by de novo sequencing using MALDI mass spectrometry. Baptide 1 comprised 7, baptides 2 and 3-10 amino acid residues, the latter being acetylated at the N-terminus. This is the first indication for the presence of such post-translational modification in snake venom proteins. None of the peptides contain cysteine residues. For biological activity studies the peptides were prepared by solid phase peptide synthesis. Baptide 3 and 2 blocked acetylcholine-elicited currents in isolated Lymnaea stagnalis neurons with IC50 of about 50 µM and 250 µM, respectively. In addition baptide 2 blocked acetylcholine-induced currents in muscle nAChR heterologously expressed in Xenopus oocytes with IC50 of about 3 µM. The peptides did not compete with radioactive α-bungarotoxin for binding to Torpedo and α7 nAChRs at concentration up to 200 µM that suggests non-competitive mode of inhibition. Calcium imaging studies on α7 and muscle nAChRs heterologously expressed in mouse neuroblastoma Neuro2a cells showed that on α7 receptor baptide 2 inhibited acetylcholine-induced increasing intracellular calcium concentration with IC50 of 20.6 ± 3.93 µM. On both α7 and muscle nAChRs the suppression of maximal response to acetylcholine by about 50% was observed at baptide 2 concentration of 25 µM, the value being close to IC50 on α7 nAChR. These data are in accord with non-competitive inhibition as follows from α-bungarotoxin binding experiments. The described peptides are the shortest peptides without disulfide bridges isolated from animal venom and capable to inhibit nAChR by non-competitive way.

Quantitative Proteomic Analysis of Venoms from Russian Vipers of Pelias Group: Phospholipases A2 are the Main Venom Components.

Venoms of most Russian viper species are poorly characterized. Here, by quantitative chromato-mass-spectrometry, we analyzed protein and peptide compositions of venoms from four Vipera species (V. kaznakovi, V. renardi, V. orlovi and V. nikolskii) inhabiting different regions of Russia. In all these species, the main components were phospholipases A2, their content ranging from 24% in V. orlovi to 65% in V. nikolskii. Altogether, enzyme content in venom of V. nikolskii reached ~85%. Among the non-enzymatic proteins, the most abundant were disintegrins (14%) in the V. renardi venom, C-type lectin like (12.5%) in V. kaznakovi, cysteine-rich venom proteins (12%) in V. orlovi and venom endothelial growth factors (8%) in V. nikolskii. In total, 210 proteins and 512 endogenous peptides were identified in the four viper venoms. They represented 14 snake venom protein families, most of which were found in the venoms of Vipera snakes previously. However, phospholipase B and nucleotide degrading enzymes were reported here for the first time. Compositions of V. kaznakovi and V. orlovi venoms were described for the first time and showed the greatest similarity among the four venoms studied, which probably reflected close relationship between these species within the "kaznakovi" complex.